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A search of neutrino magnetic moments was carried out at the Kuo-Sheng Nuclear Power Station 
at a distance of 28 m from the 2.9 GW reactor core. With a high purity germanium detector of 
mass 1.06 kg surrounded by scintillating Nal(Tl) and CsI(Tl) crystals as anti-Compton detectors, 
a detection threshold of 5 keV and a background level of 1 kg _1 keV _1 day _1 near threshold were 
achieved. Details of the reactor neutrino source, experimental hardware, background understanding 
and analysis methods are presented. Based on 570.7 and 127.8 days of Reactor ON and OFF data, 
respectively, at an average Reactor ON electron anti-neutrino flux of 6.4 x 10 12 cm~ 2 s~ 1 , the limit 
on the neutrino magnetic moments of fiu e < 7.4 x 10~ n hb at 90% confidence level was derived. 
Indirect bounds on the v e radiative decay lifetimes were inferred. 

PACS numbers: 14.60.Lm, 13.15.+g, 13.40.Em 



I. INTRODUCTION 

The strong evidence of neutrino oscillations from the 
solar, atmospheric as well as long baseline accelerator 
and reactor neutrino measurements implies finite neu- 
trino masses and mixings [3, 0]. Their physical origin 
and experimental consequences are not fully understood. 
Experimental studies on the neutrino properties and in- 
teractions can shed light to these fundamental questions 
and provide constraints to the interpretations in the fu- 
ture precision oscillation experiments. 

The couplings of neutrinos with the photons are 
generic consequences of finite neutrino masses, and are 
one of the important intrinsic neutrino properties Q 
to explore. The neutrino electromagnetic vertex can 
be parametrized by terms corresponding to interactions 
without and with its spin, identified as the "neutrino 
charge radius" and "neutrino magnetic moments" , re- 
spectively. 

This article reports on a search of the neutrino mag- 
netic moments with reactor neutrinos at the Kuo-Sheng 
(KS) Nuclear Power Station in Taiwan. R extends over 
our previous publication 0] and covers the reactor neu- 
trino spectra (4>(i7 B )) and the other experimental features 



in details. A factor of three larger data sample was used, 
and an analysis procedure combining information from 
all measured spectra was devised. 

II. NEUTRINO MAGNETIC MOMENTS 

An overview on the particle physics aspects of neutrino 
magnetic moments can be referred to a recent review [f| 
and the references therein. The neutrino magnetic mo- 
ment {nv) is an experimentally observable parameter 
which characterizes a possible coupling between neutrino 
mass eigenstates i>i and i/j with the photon whereby the 
helicity-state is flipped: 

Oi)i - 7 - {vj)r (1) 

as depicted schematically in Figure[T] The vertex marked 
"?" denotes the unknown physics to be explored. The 
parameter /i„ is usually expressed in units of the Bohr 
magneton: 

fJ-B = 7; — ; c 2 = 47ra om , (2) 
2m e 

where a em is the fine-structure constant and m e is the 
electron mass. Both diagonal and transition moments 
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are possible, corresponding to the cases where i = j 
and i ^ j, respectively. Symmetry principles place con- 
straints on the possible channels Q and require that the 
diagonal moments vanish for Majorana neutrinos. The 
study of neutrino magnetic moments is, in principle, a 
way to distinguish between Dirac and Majorana neu- 
trinos — a crucial unresolved issue in neutrino physics. 
For example, if a finite \x v is observed in the laboratory, 
and electron anti-neutrinos (z7 e ) from the Sun is detected 
with a spectrum consistent with the oscillation and solar 
model parameters, neutrinos would be Majorana parti- 
cles. Recent derivations of model-independent \i v - 
ranges indicated that upper bounds for Dirac neutrinos 
are several orders of magnitude more stringent than the 
current experimental limits. Consequently, observations 
of /i„ at the present sensitivities will imply that neutrinos 
are Majorana. 



FIG. 1: The schematic diagram of the neutrino-photon inter- 
action vertex which involves a change of the helicity-states. 
The coupling is parametrized by the neutrino magnetic mo- 
ment. 

Once a model is specified, the vertex "?" is known and 
fi v can be calculated from first principles. Minimally- 
Extended Standard Model with massive Dirac neutri- 
nos § gives \x v ~ lQ~ w \m v /l eV] /ib which is far 
too small to have any observable consequences. Incor- 
poration of additional physics, such as Majorana neu- 
trino transition moments or right-handed weak currents, 
can significantly enhance \i v to the experimentally rel- 
evant ranges [9(. Supersymmetry as well as extra- 
dimensions can also contribute to the process. 

The parameter \x v for neutrinos with energy E v pro- 
duced as vi at the source and after traversing a distance 
L is given by 11 1 



/i I/ 2 (i/ i ,L,E 1/ )=^|^U K e" 1 



(3) 



where Ujj is the neutrino mixing matrix and /iy are the 
coupling constants between v\ and Uj with the photon. 
Consequently, fx u is therefore an effective and convoluted 



parameter and the interpretations of experimental results 
depend on the exact v\ compositions at the detectors. In 
this paper, we write /i^ c = k c x 10~ 10 /ib and work with 
the parameter k c for simplicity. 

The neutrino spin-flavor precession (SFP) mecha- 
nism with or without matter resonance effects in 
the solar medium, has been used to explain solar neutrino 
deficit [l3| . This scenario is compatible with all solar neu- 
trino data. The terrestrial KamLAND experiment, how- 
ever, recently confirmed the Large Mixing Angle (LMA) 
parameter space of the matter oscillation scenario as the 
solution for the solar neutrino problem 14j . such that 
SFP can be excluded as the dominant contribution in 
solar neutrino physics fljj |. Alternatively, the measured 
solar neutrino z/ Q -e spectral shape has been used to set 
limit of Kq < 1.1 at 90% confidence level (CL) for the "ef- 



fective" vq magnetic moment (111 Il6l | which is different 
from that of a pure i7 c state derived in reactor experi- 
ments. Other astrophysical bounds on [l v were mostly 
derived from the consequences from a change of the 
neutrino spin-states in the astrophysical medium 0, 
typically of the range K as tro < 10 2 — 10 3 . However, 
these bounds are model-dependent and involve implicit 
assumptions on the neutrino properties. 

Direct laboratory experiments [5[ on /i„ utilize acceler- 
ator and reactor neutrinos as sources, and are conducted 
under controlled conditions. The most sensitive searches 
are usually performed by experiments studying neutrino- 
electron scatterings (l8| : 



v\ + e~ 



(4) 



The experimental observable is the kinetic energy of the 
recoil electrons (T). A finite fi^ will contribute to a dif- 
ferential cross-section term given by: 



da 

dT " 



2 l-T/£„ 



T 



(•5) 



The signature is an excess of events over those due to 
Standard Model of electroweak interactions (SM) and 
other background processes, which exhibit the charac- 
teristic 1/T spectral dependence. Limits from negative 
searches are valid for both Dirac and Majorana neutrinos 
and for both diagonal and transitional moments. 
The neutrino radiative decay [l9| for the process 



v s + 7 



(6) 



is another manifestation of the neutrino electromagnetic 
couplings where a change of the neutrino helicity-states 
takes place. A final-state real photon is produced in the 
process, unlike the /i„-cffccts on neutrino-electron scat- 
terings, where only virtual photons are involved. The 
decay rate and the decay lifetime is related to /iy 



r = 



i K 2 



8tt 



/4 



(7) 
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where rriij are the masses for i/jj. 



III. DIRECT SEARCHES WITH REACTOR 
NEUTRINOS 

Reactor neutrinos provide a sensitive probe for labo- 
ratory searches of fj,ir B , taking advantages of the high i7 c 
flux, low Ej/ and better experimental control via the reac- 
tor ON/OFF comparison. The measurable electron recoil 
energy spectra due to SM and [i v (at k g =1.0) scatterings 
with reactor-i7 e 's, denoted respectively by $ e (5M) and 
$ e (^„), are displayed in Figure A finite (i ffe would 
manifest itself as excess of events in the ON spectra over 
the background as derived from the OFF data, which 
have a 1 /T energy spectrum. 



was set with 66.6 days of data at a threshold of 700 keV. 
The lowest-energy bin was at a ~ 2a excess over the SM 
value. 
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FIG. 2: The expected recoil spectra due to neutrino-electron 
scattering with reactor u c at a flux of 10 13 cm -2 s -1 . Con- 
tributions from Standard Model [$ e (5M)] and magnetic mo- 
ments [<E> e (pi„)] at 10 -10 [is are shown. 

Neutrino-electron scatterings were first observed in the 
pioneering experiment [2~0 | at Savannah River where plas- 
tic scintillators were adopted as target surrounded by 
Nal(Tl) crystal scintillators as anti-coincidence counters. 
A revised analysis of the data by Ref [l8| with improved 
input parameters gave a positive signature consistent 
with the interpretation of a finite /v c at k c —2~ 4. An 
intrinsic error source of this experiment was the use of 
a proton-rich target. The cross-section for the proton- 
capture reaction: 



P 



(8) 



is much larger than that for i7 e -e scatterings. Other re- 
sults came from the Kurtchatov and Rovno [23| ex- 
periments which quoted limits of n c < 2.4 and < 1.9 at 
90% CL, respectively. A recent experiment MUNU [24 1 
at the Bugey reactor adopted a time projection cham- 
ber with CF4 gas surrounded by active liquid scintillator 
as anti-Compton vetos. A limit of k c < 0.9 at 90% CL 



A global analysis was performed 25j which combined 
simultaneously the \i v data from the reactor and solar 
neutrino experiments, as well as the LMA oscillation pa- 
rameters constrained by solar neutrino and KamLAND 
results. Only Majorana neutrinos were considered such 
that there were only transition moments. A "total" mag- 
netic moment vector A=(/Z23, M3i; M12) was defined, such 
that its amplitude was given by \A\ 2 = ^Tr(fi + fi). A 
global fit produced 90% CL limits of |A| < 4.0 x 10~ 10 fi B 
from solar and KamLAND data only, while |A| < 1.8 x 
10~ 10 /iB when reactor data were added. The results in- 
dicate the important role of reactor experiments in con- 
straining the magnetic moment effects. 

The systematic effects related to the poorly-known 
neutrino spectra were difficult to control in the 1 MeV 
range relevant to the previous experiments [26| . The ap- 
proach of the KS experiment is to achieve a threshold of 
~10 keV range through a matured and reliable detector 
technology — high-purity germanium (HPGe) detector. 
Three important advantages can be realized which sig- 
nificantly enhance the sensitivities and robustness of the 
results: (1) the potential signal rate is much increased 
due to the 1/T energy dependence of Eq.[5j (2) as shown 
in Figure O the event rates from $ e (/z„) are 24 times 
those of $ e (SM) at T~10 keV and K e ~l, such that the 
uncertainties in the SM background do not affect the \x v 
signals; (3) Eq. \5\ is mostly independent of E„ at ~10- 
100 keV, such that the ^ signals depends only on the 
well-known total reactor neutrino flux but not the details 
of 4>{i7 e ), thereby reducing the systematic uncertainties. 
In contrast, the poorly-modeled 4>(v c ) at ~1 MeV gives 
rise to uncertainties in the irreducible $ e (SM) back- 
ground (which is 2.7 times the <£> e (/i„)-signals at K e ~ 1) 
for experiments operating at this energy range. 

We report in this article data taken at KS with HPGe. 
There are three Reactor ON/OFF data taking periods, 
the key information of which are summarized in Table [U 



IV. REACTOR NEUTRINO SPECTRUM 

The i7 e 's emitted in power reactors are predominantly 
produced through /3-decays of (a) the fission products, 
following the fission of the four dominant fissile isotopes: 
235 U, 238 U, 239 Pu and 241 Pu, and (b) 239 U, following the 
neutron capture on the 238 U fuel: 238 U(n,7) 239 U . 

The KS neutrino laboratory, as depicted schematically 
in Figure [3l is located at a distance of 28 m from Core 
#1 of the Kuo-Sheng Nuclear Power Station in Taiwan. 
The nominal thermal power output is 2.9 GW. The stan- 
dard operation includes about 18 months of Reactor ON 
time at nominal power followed by about 50 days of Re- 
actor outage OFF period when about a third of the fuel 
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TABLE I: Summary of the key information on the three data taking periods. 



Period 


Data Taking 


Reactor ON 


Reactor ON 


Reactor OFF 


Reactor OFF 


DAQ 


Average 




Calendar Time 


Real Time 


Live Time 


Real Time 


Live Time 


Live Time 


i> c flux 






(days) 


(days) 


(days) 


(days) 


(96) 


(10 12 cm^s" 1 ) 


I 


July 2001 - April 2002 


188.2 


180.1 


55.1 


52.7 


95.7 


6.29 


II 


Sept. 2002 - April 2003 


125.8 


111.7 


34.4 


31.5 


89.4 


6.53 


III 


Sept. 2004 - Oct. 2005 


303.9 


278.9 


48.7 


43.6 


91.5 


6.51 



Total - 617.9 570.7 138.2 127.8 92.4 6.44 



TABLE II: Fractional mass compositions, relative fission contributions and neutrino yield per fission for the four fissile isotopes 
and the 238 U neutron capture. All values are time-dependent and only typical levels are given. 



Channels 


Fractional Compositions 
by Mass (%) 


Relative Rates 
per Fission 


Neutrino Yield 
per Event 


Neutrino Yield 
per Fission 


23b U Fission 


1.5 


0.55 


6.14 


3.4 


238 U Fission 


98.0 


0.07 


7.08 


0.5 


239 Pu Fission 


0.4 


0.32 


5.58 


1.8 


241 Pu Fission 


<0.1 


0.06 


6.42 


0.4 


238 U (n, 7 ) 239 U 




0.60 


2.00 


1.2 


Total 








7.3 



Kuo-Sheng Nuclear Power Station : Reactor Building 




FIG. 3: Schematic layout of the Kuo-Sheng Neutrino Labo- 
ratory in relation to the core of the power reactor. 



elements are replaced. Reactor operation data on the 
thermal power output and control rod status as func- 
tions of time and locations within the core were provided 
to the experiment by the Power Station. A set of soft- 
ware programs [27| was specifically developed, in associ- 
ation with the commercial SIMULATE-3 and CASMO-3 
codes [28| , both of which are extensively used in the field 
of nuclear reactor core analyses. The variations of the 
thermal power output, as well as the fission rates and 
j/ -flux of the fissile isotopes during Period-Ill are dis- 
played in Figures |3Jl-c. At steady state operation, the 
total fission rates and the total neutrino fluxes were con- 
stant to better than 0.1% and 0.2%, respectively. Data 
taken during the short durations of unscheduled reactor 
stops were included into the Reactor-OFF category. 
The typical ;7 e -yield for these channels, as well as their 



relative contributions per fission are summarized in Ta- 
blc|TTJ The ty-yields for the fission elements were adopted 
from Ref. [29[ , following a survey of the fission daughter 
isotopes and the subsequent /3-decays necessary to reach 
stability. The 238 TJ neutron capture rate of 0.6 per fission 
were evaluated by two independent methods: (a) via full 
neutron trans por t calculations of the neutrons in the re- 
actor core [30L l3ll] , and (b) by evaluating the difference in 
the decrease of the amount of 238 U and the 238 U fission 
rate. Results derived by both methods were consistent 
to a few percent. 

The <j>(v e ) 's of the five channels, as depicted in Fig- 
ure EH, were adopted from Ref. 18] for the fission ;7 c 's, 
while those following 238 U neutron capture were derived 
from standard /?-spectra of 239 U and 240 Np. The com- 
ponents were summed according to the relative contribu- 
tions per fission, and the resulting total 4>(y Q ) is shown 
in Figure Eh . This spectrum was used as input in deriv- 
ing the expected electron recoil spectra in Figure O The 
evaluated j7 c -fluxes for the three periods are given in Ta- 
ble HI where the weighted average is 6.4 x 10 12 cm _2 s -1 . 

The modeling of (j>{v c ) above the threshold of the re- 
action of Eq. [8] from reactor operation data has been 
well-established. Accuracies of better than 1.4% and 5% 
between calculations and measurements were achieved in 
the integrated [32| and differential [33| spectra. The eval- 
uation of 4>(v c ) at energy below 2 MeV, on the other 
hand, is much more complicated [26} . Many input param- 
eters remain unknown and there are no measurements to 
cross-check. Consequently, the & e (SM) recoil spectra 
below the MeV range in Figure [2] were subjected to large 
uncertainties. However, the spectra ^ e (fx v ) at 10-100 keV 
were accurately predicted due to the E„-independence of 
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FIG. 4: The variations of the (a) thermal power output, as 
well as the (b) relative fission rates and (c) i7 c -nux of the fissile 
isotopes over Period-Ill of data taking. The reactor outage 
OFF period was completed on Day-0. 



Eq.[3 They followed a 1/T profile and depended only on 
the total z7 -fluxes of Table HI which were calculated with 
an expected accuracy of a few percent. 



FIG. 5: (a) Spectral shape for reactor neutrinos due to indi- 
vidual production channels, (b) Total spectrum at the typical 
power reactor operation. 



V. EXPERIMENTAL SET-UP 

A research program on low energy neutrino physics [34j 
is being pursued by the TEXONO Collaboration at 
the KS Neutrino Laboratory. The laboratory is 
equipped with an outer 50-ton shielding structure de- 
picted schematically in Figure [5J consisting of, from 
outside in, 2.5 cm thick plastic scintillator panels with 
photo-multiplier tubes (PMTs) readout for cosmic-ray 
veto (CRV), 15 cm of lead, 5 cm of stainless steel support 
structures, 25 cm of boron-loaded polyethylene and 5 cm 
of OFHC copper. The innermost volume with a dimen- 
sion of 100x80x75 cm 3 provides the flexibilities of plac- 
ing different detectors for different physics topics. During 
the I-III data taking periods listed in Table U both the 
HPGe and a CsI(Tl) scintillating crystal array [35| to- 
gether with their associated inner shieldings were placed 
in the inner volume. The CsI(Tl) array is for the mea- 
surement of neutrino-electron scattering cross-sections. 
The /^„-search reported in this article was performed with 
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the HPGe detector. 
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FIG. 6: The shielding design of the KS Neutrino Laboratory. 
Detectors and inner shieldings were placed in the inner target 
volume. 

As indicated in Figure [31 the laboratory is located at 
the ground floor of the reactor building at a depth of 
12 m below sea-level and with about 25 m water equiv- 
alence of overburden. The primary cosmic-ray hadronic 
components are eliminated while the muon flux is re- 
duced by a factor of 4. Ambient 7-background at the 
reactor site is about 10 times higher in the MeV range 
than that of a typical laboratory, dominated by activity 
due to 60 Co and 54 Mn (half-lives 5.27 years and 312 days, 
respectively) present as dust in the environment. Both 
isotopes are produced by neutron activation on the con- 
struction materials at the reactor core. The dust can 
get settled on exposed surfaces within hours and is dif- 
ficult to remove. Some earlier prototype detectors were 
contaminated by such. Attempts to clean the surfaces 
in situ resulted in higher contaminations. Accordingly, 
the various detector and inner shielding components were 
carefully washed and wrapped by several layers of plas- 
tic sheets before transportation to KS. Outer layers were 
removed in situ only prior to installation, while inner lay- 
ers were replaced between data taking periods. Neutron 
background is comparable to that of a typical surface lo- 
cation, with no observable differences between the ON 
and OFF periods. 

The HPGe set-up is shown schematically in Figure [7] 
It is a coaxial germanium detector [37| with an active 
target mass of 1.06 kg. The construction materials and 
detector geometry followed the "Ultra Low Background 
(ULB)" standards. The lithium-diffused outer electrode 
is 0.7 mm thick. The end-cap cryostat, also 0.7 mm thick, 
is made of OFHC copper. Both of these features provide 
total suppression to ambient 7-background below 60 keV, 
such that events below this energy are either due to in- 
ternal activity or ambient MeV-range 7's via Compton 
scattering. Consequently, the background profile should 



be smooth and continuous below 60 keV down to the 
ranges where atomic effects are important. The HPGe 
was surrounded by an anti-Compton veto (ACV) detec- 
tor system made up of three components: (1) an Nal(Tl) 
"well-detector" of thickness 5 cm that fit onto the end- 
cap cryostat and is coupled to a 12 cm PMT readout 
through an additional 7 cm of CsI(Tl) as active light 
guide; (2) an Nal(Tl) "ring-detector" of thickness 5 cm 
at the joint of the cryostat; and (3) a 4 cm thick CsI(Tl) 
"base-detector" at the bottom. All of these ACV detec- 
tors were assembled within mechanical structures made 
of OFHC copper, and were read out by PMTs with low- 
activity glass. The assembly was further surrounded by 
3.7 cm of OFHC copper inner shielding. Another OFHC 
copper wall of thickness 10 cm provided additional shield- 
ings on the side of the liquid nitrogen dewar and pre- 
amplifier electronics. The inner shieldings and detectors 
were covered by a plastic bag connected to the exhaust 
line of the dewar, serving as a purge for the radioactive 
radon gas. 



Plastic Bag for Radon Purging 




Liquid nitrogen dewar 



OFBC Copper 



FIG. 7: Schematic layout of the HPGe with its anti-Compton 
detectors as well as inner shieldings and radon purge system. 

The electronics and data acquisition (DAQ) sys- 



tem [36| of the HPGe detector assembly is illustrated 
schematically in Figure [8] The HPGe pre-amplifier sig- 
nals were distributed to two spectroscopy amplifiers at 
the same 4 /xs shaping time but with different gain fac- 
tors. The amplifier signals were fed to a discriminator 
with a minimal threshold. The discriminator output pro- 
vided the on-line triggers, ensuring all the events down 
to the electronics noise edge of 5 keV were recorded. The 
amplifier output of the HPGe and the PMT signals from 
the ACV detectors were recorded by 20 MHz Flash Ana- 
log to Digital Convertor (FADC) modules with 8-bit dy- 
namic range for a duration of 10 fj,s and 25 fj,s before and 
after the trigger, respectively. The discriminator out- 
put for the various channels of all three systems (HPGe, 
ACV, CRV), as well as the timing output of the CRV 
PMTs, were also recorded. The redundancy of the read- 
out channels enhanced robustness and stability to the 
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detector performance. A random trigger (RT) was gen- 
erated by an external clock at a sampling rate of 0.1 Hz. 
The RT events provided accurate measurements of the 
pedestal levels, the DAQ dead time as well as the various 
efficiency factors. The data were read out via the VME 
bus through an VME-PCI interface [HI to a PC running 
on Linux operating system. The data were saved on hard 
disks. The DAQ system remained active for 2 ms after 
a trigger to record possible time-correlated signatures. 
The typical data taking rate for the HPGe sub-system 
was about 1 Hz. The DAQ dead time was about 4 ms 
and 2 ms per event for the HPGe and CsI(Tl) triggers, 
respectively. The system live times are listed in Table HI 
They varied among the three periods because of the dif- 
ferent trigger rates. 
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FIG. 8: Schematic layout of the electronics and data acquisi- 
tion systems of the HPGe and the associated ACV and CRV 
detectors. 



The KS Laboratory is connected with the control cen- 
ter at Academia Sinica (AS) through the telephone lines 
via MODEM. Internet connections are prohibited due to 
security reasons. During steady-state data taking, the 
experiment operated automatically without the necessity 
of human presence. The operating conditions were con- 
stantly monitored via the MODEM. Stability of key pa- 
rameters like temperatures, liquid nitrogen levels, trigger 
rates and DAQ live times were checked. The KS Nuclear 
Power Station can be reached by an hour's drive from 
AS. The laboratory was visited typically once a week, 
when manual checks and calibrations were performed. 
Hard disks were retrieved and brought back to AS where 
data were archived and copied. A total of ^125 Gbytes 
of raw data were recorded in the three data taking pe- 
riods by the HPGe system. The data to be processed 
were installed on to an external disk array storage sys- 
tem [3^] with 1.6 Tbytes total storage capacity. Relevant 
and high-level data on physical quantities were then ex- 
tracted and distributed among the users, on which data 
analysis were performed. 



VI. DATA ANALYSIS AND UNDERSTANDING 

Scatterings of z7 c -e inside the Ge target would man- 
ifest as "lone-events" uncorrelated with other detector 
systems. These events were extracted from the raw data 
through selection criteria including pulse shape analy- 
sis (PSA), anti-Compton (ACV) and cosmic-ray vetos 
(CRV). 



(a) 
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FIG. 9: Selection procedures of the recorded data : (a) pulse 
shape analysis, (b) anti-Compton selection, and (c) cosmic- 
ray veto. 
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As displayed in the correlation plot between pulse area 
and amplitude in Figure [9^, spurious background due to 
accidental and delayed "cascade" events were suppressed 
by PSA. In particular, the structure at 66 keV was from 
the decays of 73 Ge* via the emissions of two photons 
separated by a half-life of 4.6 /xs. The ACV and CRV 
cuts suppressed Compton scattering and cosmic-ray in- 
duced events, respectively. The scattered plot of HPGe 
and Nal(Tl) well-detector events is shown in Figure [9b- 
A software threshold of about 5 keV was adopted for 
the ACV detector to define the veto regions. The vari- 
ous slanted bands were due to full energy depositions of 
7-rays in the two detectors, such as those from 40 K and 
the 238 U and 232 Th series. The vertical band at 2614 keV 
originated from the coincidence of the 7-ray cascade fol- 
lowing the decays of 208 Tl. The timing correlations be- 
tween CRV and HPGe events are shown in Figure [5b. 
The PMT pulses from the CRV were fast (10 ns rise time) 
and therefore arrived before the HPGe signals which were 
processed through amplifiers with 4 /is shaping time. The 
time difference increased with lower energy because the 
trigger timing was defined by a constant-threshold dis- 
criminator. The selected regions for the lone-events are 
displayed in the respective plots in Figure [9j 

The background suppression factors as well as the sig- 
nal survival efficiencies for the three data taking peri- 
ods are summarized in Table Mil The RT events were 
uncorrelated to the other parts of the detector systems 
— similar to the neutrino- induced events. The survival 
probabilities of the RT events along the various stages of 
the analysis procedures provided accurate measurements 
of the DAQ and analysis efficiencies. The DAQ live time 
is the ratio between RT events actually recorded on disk 
to the total numbers of RT signals generated by the clock. 
The ACV and CRV selection efficiencies are the fractions 
of the RT events which survived the cuts. Only loose cuts 
of ±5(7 around the signal band were used in the PSA cut 
as shown in Figure [5b., such that its efficiency is close to 
unity. 

As illustrations, the measured spectra in Period-Ill 
before and after the ACV+CRV cuts for the Reactor 
ON data are depicted in Figure [lOl The spectra after 
the event selection criteria, also for Period-Ill, are dis- 
played in Figures [TTk and [Tib for the high and low en- 
ergy ranges, respectively. The main 7-lines were identi- 
fied and tabulated in Table HVl Most of the 7-activities 
were due to natural background from the 238 U and 232 Th 
series, as well as from 40 K and 235 U. The existence of 
activities from 235 U is particularly instructive. The low- 
energy lines implies that the source should be close to the 
target, most probably inside the OFHC cryostat, where 
238 U and 232 Th should also be present in larger amount. 
These activities were likely to originate from the front- 
end pre-amplificr components located in the vicinity of 
the HPGe target. These were expected to be one of the 
important background sources. 




After ACV+CRV Cuts 



20 40 60 80 100 120 

Energy (keV) 

FIG. 10: Measured spectra before and after the ACV+CRV 
cuts for the Reactor ON data in Period-Ill. 

Besides natural radioactivity, there were evidence of 
60 Co and 54 Mn, both being known contaminations within 
the reactor building. Decay profiles were observed for 
54 Mn consistent with its nominal decay life-time. The 
peak at 662 keV was due to 137 Cs activity in the CsI(Tl) 
ACV detector located near the HPGe. There was an 
unidentified line at 249.8 keV. The intensity was uniform 
to <4% in both ON and OFF periods but the short- 
duration rates fluctuated more than the statistical uncer- 
tainties would allow. No interpretations consistent with 
the other spectral features could be found. 

The various 7-lines also provided in situ energy cali- 
bration as well as independent measurements and consis- 
tency checks to the stability and efficiency factors com- 
plementary to the RT events. Details on background 
and detector stabilities will be presented in Section VII 
in connection with the discussions on the systematic un- 
certainties. 

The low energy spectrum of Figure lllb is relevant 
to the studies of neutrino magnetic moments. A de- 
tector threshold of 5 keV and a background of ~1 
keV _1 kg _1 day _1 above 12 keV were achieved. The back- 
ground level is comparable to the typical range in under- 
ground Cold Dark Matter experiments. This is a notable 
achievement for an experiment at shallow-depth — and 
in the vicinity of a power reactor core. Several 7-lines 
can be identified: (a) Ga X-rays at 10.37 keV and 73 Ge* 
at 66.7 keV from internal cosmic-induced activities, (b) 
234 Th at 63.3 keV and 92.6 keV from the 238 U series due 
to residual ambient radioactivity close to the target, and 
(c) a line at 80.9 keV in the OFF spectrum, to be exam- 
ined in a subsequent paragraph. The copper cryostat as 
well as the inactive surface electrode of the HPGe pro- 
vided total suppression to external low energy photons. 
Consequently, the measured spectra did not exhibit any 
structures above the the Ga X-rays end point at 12 keV 
up to about 60 keV. 

Time variations of two features were observed in the 
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TABLE III: Summary of the event selection procedures as well as their background suppression and signal efficiency factors. 
The survival probabilities are defined as the ratio between the events selected to those prior to these cuts but after previous 
cuts were applied. 



Event Selection 




Background Suppression 




Signal Efficiency 


Period 


I 


II 


III 


I 


II 


III 


Raw Data 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Pulse Shape Analysis (PSA) 


>0.99 


>0.99 


>0.99 


>0.99 


>0.99 


>0.99 


Anti-Compton Veto (ACV) 


0.054 


0.051 


0.058 


0.99 


0.99 


0.99 


Cosmic-Ray Veto (CRV) 


0.92 


0.85 


0.80 


0.95 


0.94 


0.93 


Combined Efficiency 


0.050 


0.043 


0.046 


0.95 


0.93 


0.92 



(a) 



(a) 




o10 - 



1000 1500 2000 
(b) 
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Energy (keV) 




100 120 
Energy (keV) 

FIG. 11: Measured spectra in Period-Ill after all the selec- 
tion criteria were applied for: (a) the full energy range and 
(b) the range below 120 keV relevant to magnetic moment 
studies. Key 7- lines were identified. Both Reactor ON and 
OFF spectra are separately displayed in (b). 



low energy spectrum displayed in Figure II lb . The in- 
tensity of the 10.37 keV Ga X-ray peak decayed with 
time, as illustrated in Figure [Elk in the case of data 
from Period-Ill. The best-fit half-lives of 275±14 days 
and 9.6±3.0 days to a double-exponential function agree 
well with the decays of 68 Ge and 71 Ge both of which were 
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FIG. 12: The time-variation plots in Period-Ill data taking 
for 7-lines due to (a) Ga X-rays at 10.37 keV and (b) the 
80.9 keV peak interpreted as 7's from 133 Xe. 



neut ron- act i vat ed : 
70 Ge(n,3n) 68 Ge : 68 Ge + e" 

and 



> 68 Ga (ti = 270.8 days) 
(9) 



70 Ge(n, 7 ) n Ge : 71 Ge + e" -> 71 Ga (n = 11.4 days) . 

('") 

The observed exponential decays indicated that there was 
less neutron activation on the Ge-target inside the shield- 
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TABLE IV: Summary of 7-lines intensity measured in Period- 
Ill. 



Rnerffv 


Isotopes 


Source / 


7- 1 

2 


Intensity 


(kev) 




Decay Series 




(kg day ) 


66.7 


Ge 


cosmic 


0.5 s 


15.4±0.4 


92.6 


234 mi 

Th 


238 tj 


24.1 d 


11.9±0.5 


143.8 


235 TT 
U 


235tt 
U 


7.0x 10 y 


5.1±0.8 


185.7 
186.2 


235tt 

u 

226 d 

Ra 


235tt 
U 

238 tt 
U 


7.0x 10 y 
1600 y 


} 17.2±0.4 


238.6 


Fb 


232 rj^T 

In 


10.6 h 


18.8±0.5 


0/1 s 




unidentified 




ll.OztU.O 


295.2 


Pb 


238 tt 

U 


26.8 m 


6.3±0.3 


338.3 


228 a 

Ac 


232 ml 

Th 


6.2 h 


3.7±0.5 


351.9 


214p b 


238 jj 


26.8 m 


17.1±0.4 


463.0 


228 A „ 
Ac 


232 ml 

In 


6.2 h 


1.6±0.3 


080.2 


208 rp, 
11 


232 mi. 

In 


o.l m 


1 a a _i_n 


609.3 


214 D - 
Bi 


238 TT 
U 


19.9 m 


8.1±0.2 


661.7 


Cs 


Csl(llJ 


30.1 y 


4.6±0.2 


727.3 


Bi 


232 ml 

Th 


66.6 m 


6.4±0.2 


766.4 


234m p a 


238 jj 


1.2 m 


5.0±0.3 


785.4 
786.0 


212rj- 

Bi 
Fb 


232 ml 

Ih 

238 tt 
U 


66.6 m 
26.8 m 


} 1.7±0.4 


1 9o.O 


228 a „ 
Ac 


232 rpV 

1 h 


p. u 
6.2 h 


2. ( ±U.O 


834.8 


Mn 


reactor 


312.3 d 


3.6±0.3 


860.6 


208 m, 
11 


232 rpV 

lh 


1 

0. 1 m 


O.O±0.O 


911.2 


228 a „ 
Ac 


232 rpr 

lh 


6.2 h 


19.1±0.3 


964.8 

969.0 


228 a _ 
Ac 

228 a „ 
Ac 


232 rpV 

lh 

232 rpr 

lh 


p. u 
6.2 h 

6.2 h 


} 14.4±0.3 


1001.0 


234m -ri 

Fa 


238 tt 
U 


1.2 m 


11.4±0.3 


1120.3 


214t3- 
h>i 


238 tt 
U 


19.9 m 


6.7±0.5 


1173.2 


Co 


reactor 


5.3 y 


26.0±0.3 


1238.1 


214 rj- 

Bi 


238 tt 
U 


19.9 m 


1.2±0.2 


1332.5 


Co 


reactor 


5.3 y 


27.0±0.3 


1377.7 


214 rj- 

Bi 


238 tt 
U 


19.9 m 


1.9±0.3 


1460.8 


40 T ^ 

K 


natural 


1.3x10 s y 


106.4±1.0 


1509.2 


214 rj- 

Bi 


238 tt 
U 


19.9 m 


0.6±0.1 


1588.2 


228 a 

Ac 


232 rTTi 

lh 


6.2 h 


2.5±0.1 


1620.5 


212rj- 

Bi 


232 rrri 

lh 


66.6 m 


1.6±0.1 


1630.6 


228 Ac 


232 rj-ij^ 


6.2 h 


0.6±0.1 


1729.6 


214 Bi 


238 jj 


19.9 m 


1.1±0.1 


1764.5 


214 Bi 


238 jj 


19.9 m 


5.9±0.9 


1847.4 


214 Bi 


238 tj 


19.9 m 


0.7±0.3 


2118.6 


214 Bi 


238 jj 


19.9 m 


0.2±0.1 


2204.2 


214 Bi 


238 jj 


19.9 m 


2.3±0.4 


2447.9 


214 Bi 


238 jj 


19.9 m 


0.5±0.1 


2614.5 


208 rpj 


232 rj-ij^ 


3.1 m 


14.5±0.2 



best-fit half-life of 10.6±3.3 days is consistent with the 
interpretation of 7-emissions following /3-decays of 133 Xe 
(r^=5.24 days). This isotope is a fission fragment with 
large cumulative yield (~7%). It exists in gaseous form 
and has long enough half-life to leave the reactor core. 
Air pockets contaminated with 133 Xe might have been 
trapped in the vicinity of the detector during installation. 
Photons at this low energy were fully absorbed absorbed 
in the HPGe, such the time variations were confined to 
the peak region. 

The comparisons of the key features in the measured 
spectra from the three periods are given in Table [V] 
There were no hardware intervention to the shielding 
structures and the detectors during data taking within 
one period to ensure stability of the ambient conditions 
and of the detector operations. In between the data tak- 
ing periods, there were maintenance, hardware improve- 
ment as well as detector installation efforts. The shield- 
ing door was opened and the radon purge system was 
temporarily disconnected. The passive shielding configu- 
rations in the inner target volume were re-assembled. All 
these operations were performed during Reactor ON sit- 
uations to avoid the possibilities of contaminations in re- 
actor outage. There were no incidents of sudden surge of 
any background during the ON-to-OFF transitions from 
both the HPGe detector as well as the CsI(Tl) array in 
all the three periods. 

It can be seen from Table |V] that the 40 K activity and 
the before-cut background level at 30 keV associated with 
the experimental hardware were maintained stable dur- 
ing the different data taking periods. The HPGe system 
had a hardware failure at the end of Period-II data tak- 
ing, which required the detector to be shipped back and 
forth by air flights to the manufacturer for repairs. Con- 
sequently, there was a sharp rise in the cosmic-induced 
Ga X-ray intensity between Periods II and III. The re- 
duction in the 208 Tl background between Periods II and 
III might be due to the replacement of pre-amplifier elec- 
tronic components necessary for the repairs. The back- 
ground due to the 60 Co and 54 Mn contaminations at the 
reactor were much suppressed by the shieldings and ex- 
perimental precautions, such that the detector-associated 
40 K became the dominant 7-activity. 



ing structures at KS compared to the AS laboratory 
where the detector hardware was prepared and tested 
unshielded prior to installation. Coupled with these pro- 
cesses were the (3 + decays of 68 Ga (Q=2.9 MeV with 
86% branching ratio). They were associated with two 
511 keV photons emitted in coincidence, such that the 
ACV suppression was large. 

Another interesting feature is the presence of a line 
at 80.9 keV observed only once in the beginning of the 
Period-Ill data taking, which was a Reactor OFF dura- 
tion. The time evolution is displayed in Figure [T2"b . The 



VII. NEUTRINO MAGNETIC MOMENTS AND 
RADIATIVE DECAY SEARCHES 

Only the low energy spectra like the one shown in Fig- 
ure [TTJd were used in the studies of Since data taking 
conditions were not identical among the three periods, 
the data were analyzed and Reactor ON/OFF spectra 
were compared independently within each period. Con- 
sequently, the variations in the key background features 
shown in Table|V]do not affect the evaluations of fi„. The 
analysis procedures were extended from those used previ- 
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TABLE V: Summary of the major spectral measurements in the three data taking periods. The intensities of the Ga X-rays 
are those at the start of the data taking. The others are time-averaged over an entire period. 



Period 


Ga X-Rays 


Event Rates (kg May" 




54 Mn 


30 keV Band (kg~ 


1 keV~ 1 day -1 ) 


40 K 


208 rpj 


60 Co 


Before-Cut 


After-Cut 


I 


134±6 


107.2±0.7 


29.2±0.3 


9.8±1.0 


4.2+0.7 


31.6±0.1 


1.58±0.02 


II 


45±3 


106.4±0.9 


25.6±0.4 


14.5±1.7 


3.1±0.4 


31.2±0.1 


1.34±0.02 


III 


128±7 


106.4±1.0 


14.5±0.2 


26.5±0.2 


3.6±0.3 


30.9±0.1 


1.42±0.01 


RMS/Mean 


0.40 


4e-3 


0.27 


0.42 


0.12 


9e-3 


0.07 



ously [4J through a global treatment which incorporated 
additional constraints and measurements. 

The PSA cut identified events with spurious and convo- 
luted pulse shapes, some of which were due to electronic 
noise. These were not stable with time in general. The 
ACV and CRV cuts, however, rejected background due 
to physical processes such that the suppression factors 
should be stable during data taking provided that the 
hardware were properly controlled and cross-monitored. 
We denote the spectra surviving the PSA cut by 
and <j) l + for the Reactor OFF and ON data, respectively, 
while those after further ACV+CRV cuts similarly by 
(f>L and 4>+- The before and after-cut Reactor ON spec- 
tra </>\ and 4> + in Period-Ill are displayed in Figure fTOl 
The neutrino-induced contributions to the electron recoil 
spectra can be described by 

v = * e (5M) + k c 2 ■ Q! e (fi u = KT 10 mb) ■ (11) 

The K e 2 dependence on the /^-signals follows from Eq. [5] 
There are two constraints relating the four measured 
spectra within one period — 

(I) the excesses of 4>+ over <fr_, if any, are neutrino- 

induced: 

(j) f + = e-v + cj)L (12) 

where e is the selection efficiency of the ACV+CRV 
cuts as given in Table HTfl and 

(II) the suppression factors for the ACV+CRV cuts 
should be constant between the Reactor-ON and 
OFF data taking: 



Global minimum-^ 2 analyses were performed to the 
data from the three periods independently. Though the 
DAQ threshold was 5 keV, the analysis threshold of 
12 keV was chosen to avoid complications of the atomic 
effects 40] as well as the time- varying background from 
Ga X-rays. The background profile was continuous and 
stable from this energy threshold up to 61 keV as ex- 
pected, and could be described by a fit of a polynomial 
function to <jr_. The best- fit values of k c 2 as well as the 



X 2 /dof for the three periods are listed in Table IVT1 It can 
be seen that Eq s . [T2l and IT51 provide excellent descriptions 
to the data. 



TABLE VI: Summary of the best-fit results for the three pe- 
riods and the xVdof indicating the goodness-of-fit. 



Period 


Best-Fit 

2 


X 2 /dof 


I 


-0.52+1.05 


90/97 


II 


0.06+1.21 


108/97 


III 


-0.84+0.87 


84/97 


Combined 


-0.53+0.59 





The data taking durations of the three periods lasted 
from eight to twelve months. Demonstrations of the 
background and detector stabilities are therefore cru- 
cial. Since each period was treated independently with 
its own Reactor ON and OFF data compared, the varia- 
tions among the different periods shown in Table [V] did 
not lead to systematic effects on the analysis. Never- 
theless, the stability (an RMS spread of 0.9%) of the 
before-cut levels at 30 keV among the different periods 
suggested that the overall background was in good con- 
trol throughout the experiment. 

Two complementary studies were performed on the 
systematic effects of the experiment, and they are dis- 
cussed in details in the following paragraphs. 

TABLE VII: The measured stability levels of various key 
spectral features, based on P-III data: R± is the "resid- 
ual/statistical error" ratio for the ON/OFF data, while 5 to t 
is the error/mean ratio for the entire DAQ period. 



Monitors 


R± 


Stot 


X 2 /dof 


4U K 1462 keV Band 


-0.42 


3.5e-3 


55/53 


208 T1 2614 keV Band 


0.18 


1.3e-2 


50/53 


Band at 30 keV 


0.005 


8.2e-3 


55/53 


Band at 400 keV 


1.3 


4.1e-3 


79/73 


Suppression factor 








(ACV & CRV) at 30 keV 


-0.21 


9.6e-3 


72/73 



(1) Stability Measurements: 

Several important spectral features were continuously 
measured to monitor and demonstrate the stabilities of 
the detector performance and background conditions. 
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These stabilities are summarized in Table IVIII for the 
350 days of Period-Ill data. The monitored features in- 
cluded the (ACV+CRV) suppression factors, as well as 
the event rates of the 40 K and 208 Tl bands, together with 
the two sampling bands at 30 keV and 400 keV. 

To quantify the stabilities, the time-evolution over the 
entire period was fit to a constant. The stability parame- 
ter Stot corresponds to the ratio of error/mean of the fits. 
The excellent values on x 2 /dof indicate that the stability 
hypothesis is valid and that the fluctuations of individual 
data points were statistical in nature. Accordingly, <5 tot 
can be taken as the stability levels of the various features 
being monitored. The parameter R± denotes the ratio of 
residuals of Reactor ON— OFF data to their respective la 
statistical accuracies. The measured results of \R±\ ~ 1 
indicate that the experimental conditions were stable be- 
tween the Reactor ON/OFF data set, such that their 
differences were consistent with statistical fluctuations, 
rather than due to systematic variations in hardware or 
background. 

(2) Bounds on Possible Instabilities: 

The various channels which may contribute to possible 
instabilities to the ON/OFF data at the low energy sig- 
nal region were identified and their effects were derived 
from actual data and complete detector simulations. The 
results are summarized in Table IVlTlk -c. Their effects to 
the ON/OFF comparisons at 30 keV, relative to the over- 
all after-cut background 0(cut) ~ 1.5 kg _1 keV _1 day _1 , 
are denoted by S±. 

The bounds on the ON/OFF stabilities due to event 
rate variations are shown in Table IVIIIa . The DAQ live 
times were accurately measured to 5± <5e-4 through the 
ratios of generated to recorded RT events. Within each 
period, the HPGe detector was continuously maintained 
at liquid nitrogen temperature to minimize the diffusion 
of the inactive layer which is about 1% of the active de- 
tector mass. A change of thickness of <1% implies a 
stability level of S± <le-4. 

The various background sources which may produce 
ON/OFF instabilities are listed in Table Ivrilb , where 
their fractional variations are denoted by 5(bkg). The 
evaporation rate of the liquid nitrogen in the HPGe de- 
war was measured to be stable to <1%. There were no 
observable changes [<5(bkg)<5%] in the key spectral fea- 
tures in a test data set taken when this radon purge sys- 
tem was disconnected. Therefore, residual variations on 
the background due to fluctuations in radon diffusion are 
expected to be at the level of 8± <5e-4. 

There were five specific lines ( 54 Mn, 60 Co, 68 Ga, 133 Xe 
only for Period-Ill, and the 249 keV line) where varia- 
tions with time were observed or should be examined. 
Their effects on the signal region were derived from full 
simulation studies using the measured peak intensities 
0(peak) in Table [TV] and <5(bkg) between the ON/OFF 
periods as input and normalizations. The ON/OFF sta- 
bilities at the signal region are the S± entries in Ta- 



TABLE VIII: Identified candidate channels which may con- 
tribute to Reactor ON/OFF instabilities. Their impact on 
Reactor ON/OFF comparisons at the 30 keV signal region 
(S±) were evaluated from P-III data and simulations. The 
three categories are listed separately: (a) event rate changes, 
(b) background variations, and (c) instabilities in detector 
performance and analysis efficiencies. The intrinsic stabilities 
of the background channels in Table (b) are given by <5(bkg) 
while !)(MC) denotes the simulated background level at the 
signal region with keV -1 bin-width, normalized to one ob- 
served event at the peak. Entries in Table (c) include the 
stabilities of the detector systems <5(Detectors) and veto effi- 
ciencies <5(Veto), while the background fraction for each mode 
is given by fs. The overall uncertainty limit combined in 
quadratures is given in (d). 

(a) Event Rate Changes 

Channels ON/OFF Stability (S±) 

DAQ live time <5e-4 
HPGe target mass <le-4 



(b) Background Variations 



Background ON/OFF Stability Levels 



Source 


5(bkg) 


r)(MC) 


8± 


Radon diffusion 


<0.05 




<5e-4 


54 Mn 


<0.3 


2e-3 


<2e-3 


60 Co 


<0.06 


2e-3 


<4e-3 


68 Ga 


<0.4 


3e-8 


<le-6 


133 Xe (only Period-Ill) 


1 


4e-4 


<3e-3 


Line at 249 keV 


<0.04 


le-3 


<3e-4 



(c) Detector Performance and Dat 


a Analysis 


Detector Systems f B ON/OFF Stability Levels 
& Cuts ^(Detector) 5 (Veto) 5± 


Both CRV & ACV 0.64 <5e-3 
ACV only 0.30 <0.05 
CRV only 0.01 <0.1 
Surviving CRV+ACV 0.05 <0.15 


<le-6 <le-5 
<7e-4 <4e-3 
<0.01 <2e-3 
<5e-4 


(d) Combined Limit 


Combined (a) + (b) + (c) 


S± < 7e-3 



blc IVIIIb . derived from: 

6± ~ [ 0(peak) • (5(bkg) • ^(MC) ] / 0(cut) (14) 

where ^(MC) is the simulated background level at the 
low-energy signal region with keV -1 bin- width, normal- 
ized to one observed event at the peak. From physical 
expectations, the background sources were taken to be 
uniformly distributed among the space in the vicinity of 
the HPGe detector, and their emissions were isotropic. 
The HPGe target was hermetically surrounded by the 
ACV and the detector response was mostly uniform in 
all directions. Consequently, the effects due to the exact 
source locations were small and the background levels 
only varied by ~10% (RMS) if localized point sources 
were used instead. 
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The tight limits can be understood qualitatively as fol- 
lows. The weak intensity of the 54 Mn line and the long 
decay lifetime of the 60 Co lines suppressed the ON/OFF 
variations to 8± <2e-3 and 5± <4e-3, respectively. The 
decays of 68 Ga are mostly via the emissions of /3 + , such 
that the ACR was efficient as vetos in catching the two 
511 keV photons. The 80.9 keV 133 Xe line appeared only 
in Period-Ill OFF period. At this low energy, the photo- 
electric cross-section is large compared to the Compton 
effects. The residual Compton-scattered photons were 
of lower energies and had large absorption by the copper 
cryostat as well as by the inactive layer of the HPGe. Ac- 
cordingly, the effects of this line were confined to the peak 
region. The correlated background at the low energy 
signal regions were suppressed to the level of 5± <3e-3. 
The short-time-scale fluctuations of the unidentified line 
at 249.8 keV were washed out to better than 4% when 
the extended ON/OFF periods were integrated, giving 
5± <3e-4. 

The stabilities of the ACV and CRV detector sys- 
tems, denoted by ^(Detector), as well as their effects 
on the veto efficiencies [5 (Veto)] and 6± are summa- 
rized in Table [VTTTb . From Table [V] the raw background 
rates before the ACV+CRV cuts were about 0(raw) ~ 
30 kg _1 keV _1 day _1 at 30 keV while the after-cut back- 
ground levels were 0(cut) ~ 1.5 kg _1 keV — 1 day — 1 , such 
that the combined suppression factors were about 5%, 
as shown in Table IIIII The fractions of the raw back- 
ground suppressed by each category of selection cuts are 
denoted by fs- The ON/OFF variations at low-energy 
are therefore given by: 

S± ~ [ is ■ *(Veto) • 0(raw) ] / 0(cut) . (15) 

About 64% of the raw background were vetoed 
by BOTE ACV and CRV. These were cosmic-ray 
events which penetrated the active ACV and emitted 
bremsstrahlung photons which interacted with the HPGc 
via Compton scatterings. Besides having the redundancy 
of signals in two detector systems, the energy depositions 
of these events were of the order of tens of MeV at the 
ACV detectors, where the threshold was about 5 keV. 
Consequently, the rejection efficiencies were 100% and 
small variations in the ACV gains or thresholds had neg- 
ligible effects on the stabilities. About 30% of the raw 
background were tagged ONLY by the ACV detectors. 
These events were due to ambient radioactivity and typ- 
ically deposited 100 keV to 1 MeV energy at the ACV. 
Studies were performed with both simulated and real 
data on the effects of the changes in software threshold. 
They indicated that an instability in the ACV detector 
gain of (5(Detector)~5% would lead only to a variation 
in the rejection efficiencies of <5(Veto)<7e-4, which trans- 
lated to 6± <4e-3. Similarly, only about 1% of the raw 
background produced CRV tag alone. These were high- 
energy bremsstrahlung events induced by cosmic-ray in- 
teractions with the passive shielding materials, where the 



photons penetrated the ACV system. Since the PMTs of 
the CRV system were optimized at the stable plateau 
region, a 10% change in the PMT gain did not give 
rise to observable change in tagging efficiencies (that is, 
<J(Veto)<l%), so that 6± <2e-3. Finally, about 5% of 
the raw background survived both ACV+CRV cuts. The 
time variations of their survival probabilities were accu- 
rately measured to S± <5e-4 by studying the response to 
the same software procedures with the RT events. 



TABLE IX: The various sources of the systematic uncer- 
tainties, their magnitudes ^(Source) and their contributions 
to <5(K e 2 ). In particular, the first entry followed from Ta- 
bles [VIlMvTTl] and is the dominant contribution. 



Sources 


S (Source) 


8(k c 2 ) 


ON/OFF instabilities [S± (combined)] 


<1% 


<0.18 


Efficiencies for neutrino events 


<5e-4 


<5e-4 


Rates for $ e (ji v ) 


<3% 


<0.03 


$ e (SM) background subtraction 


23% 


0.03 


Combined Systematic Error 




<0.18 



Taking into account the measured stability levels in 
Table IVTIl at the 30 keV signal region (S to t ~9e-3), as 
well as the combined upper bounds on the possible Re- 
actor ON/OFF variations (6± <7e-3) in Table ELL^d) , 
an uncertainty of S± (combined) < 1% was adopted to 
quantify the possible Reactor ON/OFF systematic in- 
stabilities at the low energy signal region. The influence 
of this and other sources of systematic effects are sum- 
marized in Table ILX| where the stabilities of the sources 
are given by (5(Source) while their effects on the magnetic 
moment searches are represented by 5(k c 2 ). The stabili- 
ties of the efficiencies of possible neutrino-induced signals 
in Table ITTT1 were accurately measured to (5(Source)<5e-4 
through the survival probabilities of the RT events over 
the various selection criteria. The $ e (/i^) rates at 10- 
100 keV depend mostly on the total v e flux which is 
well-modeled to <2% at steady-state operation (3^|, as 
discussed in details in Section IV. Residual uncertain- 
ties arise from the finite rise and fall time of the neu- 
trino flux relative to the sharp Reactor ON/OFF 
instants, introduced through several long-lived isotopes 
(like 93 Y, 97 Zr, 132 Te) where the total fission yield is 
about 15%. About 10% of the data taking time may 
be subjected to this effect, such that a combined uncer- 
tainty of 5(Source)<3% was derived for the evaluation 
<b e (n v ). A conservative estimate of 30% systematic un- 
certainties for the <p(v e ) spectrum below 2 MeV translates 
to 23% uncertainties in the & e (SM) background rates be- 
low 100 keV. However, this only have little impact to the 
accuracies in n e 2 since the effect being studied in $ e (/i„) 
is an order of magnitude larger than $ e (SM) at 10 keV. 

Combining statistically the best-fit values on « e 2 for 
the three periods and adopting the systematic uncertain- 
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ties derived in Table llXl the results 

Kc 2 = -0.53 ± 0.59(stat.) ± 0.18(sys.) (16) 

were obtained. Adopting the unified approach [lj, a di- 
rect limit on the v e magnetic moment 



AH < 



7.4 x 1CT 11 /is 



(17) 



at 90% CL was derived. The combined residual spectrum 
for the three periods of the Reactor ON data over the 
background profiles is depicted in Figure [T3j The best- 
fit 2(T region for n e 2 is superimposed. It has been verified 
that the correct positive signals could be reconstructed 
by the same analysis procedures operating on simulated 
spectra with /z^-induced events convoluted with similar 
background profiles. The results are also not sensitive 
to alternative choices of 6± (combined) in Table IIXI given 
the constraints of the values in Tables IVIIKjVIiTl 
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FIG. 14: Summary of the results in the searches of neutrino 
magnetic moments with reactor neutrinos. Both the limits 
and the detection thresholds of the various experiments are 
shown. 
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FIG. 13: The residual plot on the Reactor ON data of all 
periods combined over the background spectra. 

Depicted in Figure [TJ] is the summary of the results 
in fi^ e searches versus the achieved threshold in reactor 
experiments. The dotted lines denote the ratios between 
^-induced and SM cross-sections [R = a(n)/a(SM)] as 
functions of (T, fiu e ). The KS(Ge) experiment has a much 
lower threshold of 12 keV compared to the other measure- 
ments. The large R-values imply that the KS results are 
robust against the uncertainties in the SM cross-sections. 
In particular, in the case where the excess events reported 
in Ref. [li| are due to unaccounted sources of neutrinos, 
the limits remain valid. 

Results from oscillation experiments [3, Q indi- 
cate that v c is predominantly a linear combination 
of mass eigenstates v\ and v-i with mixing angle #12 
given by sin 2 6*i2 ~ 0.27. The mass differences be- 
tween the mass eigenstates are Am 2 2 ~ 8 x 10~ 5 eV 2 
and Am 2 ^ ~2x 10~ 3 eV 2 . Both "normal" (nor: 
ni3 3> ni2 > mi) and "inverted" (inv: m 2 > mi 3> 1113) 
mass hierarchies are allowed. The v\ —> v% and v-i — > 1/3 



radiative decays are allowed only in the inverted mass 
hierarchy, while vi — > v\ is possible in both hierarchies. 
Adopting these as input, the limit of Eq. [XT] can be 
translated via Eq. [7] to indirect bounds of 



T13 

rrio 



{inv : v x -> V3) > 3.2 x 10 27 s/eV A (18) 



(inv : vi -> v z ) > 1.2 x 10 27 s/eV 5 (19) 
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(nor + inv : v 2 -> v\) > 5.0 x 10 31 s/eV 5 (20) 



at 90% CL. These limits are sensitive to the bare 
neutrino-photon couplings and are therefore valid for 
neutrino radiative decays in vacuum. They are summa- 
rized in Figure [15] Superimposed are the limits from the 
previous direct searches of excess 7's from reactor 41 1 
and supernova SN1987a [iH neutrinos. It can be seen 
that bounds inferred from v-e scatterings are much more 
stringent than those of the direct approaches. 



VIII. CONCLUSION AND PROSPECTS 

This article describes a study of possible neutrino- 
photon couplings using neutrinos from nuclear power re- 
actor as source. A germanium detector of target mass 
1.06 kg was adopted as target where data taking and 
analysis thresholds of 5 keV and 12 keV, respectively, 
were achieved. Sensitive limits were derived on neu- 
trino magnetic moments and neutrino radiative decay 
lifetimes. Good background level comparable to those of 
the underground dark matter experiments was achieved 
at the 10 keV range. This is the first experiment operat- 
ing with such low threshold in a high neutrino flux. The 
unique data collected with the HPGe detector at KS open 
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FIG. 15: Summary of the excluded parameter space on neu- 
trino radiative decay lifetimes from i> e measurements. Bounds 
on v\ and U2 from reactor and supernova SN1987a data are 
denoted by R and SN, respectively. The superscripts D and I 
correspond to direct limits and indirect bounds inferred from 
magnetic moment limits, respectively. The subscript "12" is 
attributed to decays driven by Am5 2 , and so on. The upper 
bound (m up ) on m„ is due to limits from direct mass mea- 
surements, while the lower bounds mf™ and mj° v are valid 
for the normal and inverted hierarchies, respectively. Bounds 
for Rj 3 and R23 are represented by the same bands in this 
scale. 



the possibilities of studying other phenomena like the in- 
trinsic properties of reactor electron neutrinos [3l| , pos- 
sible neutrino-induced nuclear transitions, and searches 
for reactor axions, all of which are being pursued. 

The sensitivities for direct searches of neutrino mag- 
netic moments [j| with neutrino-electron scatterings ex- 
periments scale as 



1 



L M t J 



(21) 



where N v is the signal events at some reference magnetic 
moments, B, M, t are the background level, detector mass 
and measurement time, respectively. The best strategy 
to enhance the sensitivities is to have an increase on N v , 
which is proportional to the neutrino flux </>„ and is re- 
lated to the detection threshold. The atomic energy level 
effects [4(j| | . however, limit the potential sensitivity im- 
provement as the threshold is reduced below the typical 
atomic scale of ~1 keV. There is an on-going experiment 
GEMMA [HI and a proposal MAMONT Q pursuing 
along these directions. 

However, this approach cannot improve on the \x v sen- 
sitivities indefinitely. Inferring from Figure [21 for reactor 
neutrinos, $ e (/•*!/) is only 2% of $ e (SM) at 1 keV as 
jtz„ — > 10~ 12 /ib- Accordingly, it would be experimen- 
tally difficult to study zz„ by looking for anomalous ef- 
fects in neutrino-electron scatterings over SM behaviour 
for [i v < 10~ 12 hb- This limitation can be evaded, at 
least conceptually, by doing the analog of "appearance" 



experiments in the case of Majorana neutrinos. One can 
look for signatures of anti-neutrinos of a different flavor in 
a pure and intense neutrino beam which passes through a 
dense medium or an strong magnetic field. Though there 
is no fundamental constraint on the lowest reach of the 
detectable vjv ratio, realistic accelerator-based experi- 
ments are still many orders of magnitude less sensitive 
than the neutrino-electron scattering experiments with 
reactor neutrinos (H)]. 

The TEXONO Collaboration meanwhile is pursuing 
an R&D program on the "ultra-low-energy" germanium 
detector. A threshold of 100 eV has been achieved with 
a 10 g prototype [4a ]. The goals are to develop a ~1 kg 
detector to perform the first experimental observation of 
neutrino-nucleus coherent scattering using reactor neu- 
trinos. The by-products of such an experiment will be 
dark matter searches at the low WIMP mass regions, as 
well as the probing of fi„ e down to the ~ 10~ n tie range. 
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